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THE GENERAL SOLUTION TO THE CLASSICAL PROBLEM OF FINITE 

EULER BERNOULLI BEAM 

M. Y. Husaalni and C. L. Amba-Rao* 

NASA Ames Research Center, Moffett Field, California 


SUMMARY 

An analytical solution is obtained for the problem of free and 
forced vibrations of a finite Euler Bernoulli beam with arbitrary 
(partially fixed) boundary conditions. The effects of linear viscous 
damping, Winkler foundation, constant axial tension, a concentrated mass, 
and an arbitrary forcing function are included in the analysis. Mo 
restriction Is placed on the values of the parameters Involved, and the 
solution presented here contains all cited previous solutions as special 
cases. 

1. INTRODUCTION 

Since the classical theory of beam was evolved by Euler and 
Bernoulli, a large literature has accumulated on the subject, the 
major part of which deals with either Infinite beams or finite beams 
with standard boundary conditions. The general solution for the dynamic 
response of the Infinite Euler-Bernoulll beam with arbitrary Initial 
conditions, subjected to an arbitrary load Including the effects of 
damping, an elastic foundation and constant axial load, was obtained 
comparatively recently by Stadler and Shreeves (1). The analogous 
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problem for a finite beam does not seem to have been solved In all 
generality. The first study allowing an elastically restrained edge 
condition, which Is of particular Interest In the aerospace field, appears 
to be that of Carmichael (2). Among subsequent Investigations, those 
of Kaltbeek (3); Hess (4), Amba-Rao (5), Sharp and Cobble (6) and 
Amba-Rao and Hussainl (7) are perhaps worth mentioning. 

Maltbaek (3) solved the problem of a uniform beam with a rigid 
mass attached at Its center and with elastic supports. Hess (4) gave 
a more complete solution. An^a-Rao (5) considered the free vibrations 
of an elastically restrained beam carrying an arbitrary nusiber of discrete 
masses. He used the transform method, treating the added masses as impulses 
(represented by Dirac delta functions) In the mass density function. 

Sharp and Cobble (6) solved the problem of the arbitrarily loaded daaq>ed 
beam elastically restrained against rotation. Amba-Rao and Hussainl (7) 
presented a closed form solution to the problem of free vibrations of a 
partially fixed, linearly damped ^ beam on a Winkler foundation, carrying 
arbitrary masses. The present note extends the latter solution to Include 
the case of forced vibrations under arbitrary initial conditions, and 
unifies all earlier investigations Into a general theory. 

2. STATEMENT OF THE PROBLEM 

In addition to the usual assumptions of the Euler-Bernoulll theory, 
the beam Is assumed to have constant axial tensile load, with viscous 
damping proportional t:o velocity, and to be resting on a Wltikler foundation 
whose resistance Is proportional to the transverse displacement. The 
relevant differential equation In operator form Is 
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JCy(x,t) - q(x,t) 


( 2 . 1 ) 


r-EJ-~-P^+f + d-^+p-^ 
ax'* 3x2 3t at* 

H«re y(x,t) Is the lateral deflection » EJ the flexural rigidity 
of the bean* P the mass density of the beam and concentrated aasses, 
X the space coordinate measured along the length of the beam, . 

P the tensile axial force on the beam, f the Winkler foundation 
constant, d the damping constant, t the time, and q the forcing 
function. The initial and boundary conditions are, respectively: 


y - u(x) , 


at t ■ 0 


1 _ 


- -R 

3x2 3x 


at X ■ 0 for all t 


1 


at X L for all t 


3x2 2 3X J 

where L Is the length of the beaijl and T^^T 2 , and R 2 are spring 
constsnts (partial fixity factors^. ' These boundary conditions are 


unmixed. The method of solution (jlVen her^ Is applicable even if the 


conditions are mixed. 


3. SOLUTION 


The general theory of linear dgwystors suggests two fundamental 


methods of solving equation (2.1). One method is to find the Inverse 




of th« oporator X, that is, in offset to obtain tha Groan's function. 

Tha oparator X is not self-adjoint with raspact to either x or t. 

However, if the Green's function G is evaluated satisfying the hoBogeneotis 
initial and boundary conditions, then the problem can be solved in 
principle by an extended definition of X as given in equation (3.1) 
below: 

Xy ■ q + (pu(x)d'(t) + ov(x)6(t) + d6(t)] 

■¥ terms due to inhomogeneous boundary conditions 
“ Q (say) (3.1) 

where £(t) is the delta function. The solution takes the form 

y(x,t) • I r G(x,t; C,x)Q(C,T)dC dt (3.2) 

In the present paper, the solution is obtained by the second method 
using spectral representation of X in terms of x-elgenfunctions. 

First, the solution Is sought as a sum of two functions, one of which is 
the solution of a homogeneous problem with Inhomogeneous boundary conditions, 
while the other is the solution of an inhomogeneous equation with homo- 
geneous boundary conditions. Thus, let 

y(x,t) - yj(x,t) -»• y^(x,t) (3.3) 


where 


5 


with the boundary conditions 

y;"(0,t) - -Tjyj(O.t) 
y;"(?l,t) - T^y/JI.t) 
y'j(o,t) - -RjyJ(O.t) 
y’j'(75,t) - R 2 y;(R,t) 

and 

^yg “ q (3.6) 

with the boundary conditions 

ygCo.t) - y2(^.t) - 0 
y'jCO.t) - y^'(3l,t) - 0 

and the initial conditions 

y^(x,0) ■ u(x) “ yj(x,0) - U(x) 

9y2 9yi 

(x,0) = v(x) - (x,0) ■ V(x) 

where, by a proper choice of the origin and unit of length, the end 
points of the beam are at x ■ 0 and x * ][, and primes denote 
differentiation with respect to x. 

a) EVALUATION OF yj 

The solution yj of equation (3.4) is given in reference 7. A 
slightly different approach is presented here which yields results 
Identical to those of reference 7. The approach is expected to be 
useful particularly if the series representing the Inverse of the finite 
Fourier transform is not summable in closed form’. 
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Write 

Vj ■ <Kx)exp(-at)co8(wt) (3.9) 

Substituting this in equation (3.4), we find that tC(x) satisfies the 
ordinary differential equation 

EJ - P ^-4 + fi|» - p(a^ + w^)^ - 0 (3.10) 

dx^ 3x^ 

where 2ap ■ d, and with concentrated masses being represented by delta 
function p “ ffl -f Mjd(x - x^). Thus, the concentrated mass is supposed to be 
situated at the point x ■ and d is assumed proportional to p. 

The boundary conditions are 

^(3)(0) - -Tiiii(O) , - Tjil-CTO 

(|i^^)(0) ■ -Rjip^^^(O) and ij»^^^(A) ■ R29^^^(A) 

where denotes 1th derivative of ip with respect to x. 

Define t 

- 1 - f . 

* 'P ■ I — f cos nx dx 
Jo 3x^ 

Integrating by parts, we establish the recurrence relation 
4^5 - ^ [(-D^ip^^'^^^K) - ip^^'^^^0)] 

- 4 - ip^^'^^^0)] + (3.12) 

n n^ 

?;(n) - a°5; - 4 [(-i)"ip^^^ 3^) - ip^^^o)] 

-4 [(-1)%^^^^) - 'P^^^o)] +4 ^'*5' (3-13) 

n^ n^ 

The expression for 4‘*lp is obtained from the finite cosine transform of 
equation (3.4): 


I 


(3.11) 
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‘"♦[i + h 

where 

Finally, 

i(n) ' 

where 

Then 

or 

4»(x) ■ 


♦(") [s + “*> - n] 


EJ J 


+ -I- S, 


+ 0)2 




^ Rn2 + 2£2){(-l)"«|»^^^a) - ♦^*^0)) 

(n2 + «,2)2 _ fc4 - 

- {(-l)"ij;^^^a) - 4»^^^(0)} + Sj] 




m 


(g2 ^ 0)2^ 

EJ 


i|»(x) - ij»(0) + ^ i(n)co8 nx 

i — r^^^M.,i.(x.) - - /^^o)) 

S £4 _ L J J 

+ Ijdx - XjD] - I [/'^0){lj(x) + (t2 + k2)l3(x)} 
+ (ijCx) + + k2)I^(x)} + {il<^^^(0)l3(x) - 


0.14) 


(3.15) 


) 


«(3)(H)IJx))] 


(3.16) 



8 


where 


Ij(x) 

Ij(x) 


1 ^ 71 cosh(n -• x)/t^ 

’ 2(&2 + k2) 2^^ + k® 8lnh 3l/l2 + i? 

1 _ X r cosh(X - x)/t^ k* . 

2(1** - k**) " 4k^ \Jl^ + k^ Blnh + k^ 


coiihCTt - x)/t^ - k^ "1 
- k^ slnh 75/i^ - k^ J 


IjCx) - Ij(7l - x) and I^(x) - 1^(75 - x) 


b) SOLUTION OF £y^ - q(n,t) 

The operator is a sum of two commutative operators and 

defined as 


L •EJ-~-P-^+f 
X ax'* 3x2 

S ■ ^ Tt 


(3.17) 

(3.18) 


The inverse of may be obtained by considering as a constant. 
The results will be a function of the operator and should be interpreted 

by using the spectral representation of L^. Consider L^: Its dooaln 

is the set of all functions in S (which is a linear vector' space of 
real-valued square Integrable functions over (0,70) and they have piecewise 
continuous fourth derivatives for 0 < x < x^ and for x^ < x < 71, whic'i 
satisfy the conditions 


'I'(O) ■ i|((7() ■ \)^"(0) ■ ij/"(7() ■ 0 


(3.19) 


and are such that belongs to S. 

If are in D^, it can be easily shown following standard 


procfcluri' 
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where the inner product le defined ae 

(*‘•*’ 4 ) ■/„ 

In other worde» is self-adjoint and ita eigenfunctions ere 
These cen be obtained from the analyaie of section (a) by letting 

Tj, Tj • and R^, R 2 "*■ 0* The defined by the relation 


(3.20) 


(3.21) 


JJT17T 


form the orthonomal basis for S. Hence q, U, V can be represented 
as 


(3.22) 


where 


q/p - Eq.(t)#.(x) 

1 ^ ^ 

OB 

U - E U^^^(x) 


V " E Vi<*) 
1 


o » ♦!<'>« 


ii ■/ «i(e.o 

•'o p 

Uj-J U(0^^(C)dt 

^ Q 


'A 


V(0*,(5)dC 


(3.23) 


(3.24) 


The equation (3.6) is written as 


» y2(*»0) ■ *’(x) , (x.O) - V(x) 


(3.25) 


I 


Trtating «i a constant, we find the solution of (3.25) (by the 
method of Laplace tranaforms or by variation of parameters) aa 

y^ - P exp[-o(t - T)] ^ .ailjtiMx 


+ exp (-at) U(«) 


+ exp(- ot) — V(x) 


(3.26) 


where S ■ - L^/p and 2 op ■ d. • It Is to be noted that 

the function of the operator la written before the function on which 

It acts and it Is Interpreted accordingly. For Instance, 


f(Lj^)U(x) 


£f(X )U ♦ < 
“ n n n 


where A are the eigenvalues of L end 6 M the associated normalised 
n X n 

eigenfunctions. Thus, 


(t) + F^(t)(oU + V ) 
i n n n 


(3.27) 


sinh 6 

F„(t) - exp(-ot) r— S- 


g and 6 ■ - X /p 

n n 
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DISCUSSION 

Substitution of (3.16) In (3.11) yields four equations for five 
unknowns 'I'^^^(O), <|i^^^(o), i|»^^^(7) and Evaluation of 

at X “ Xj from (3.16) provides the necessary fifth equation. The determinant 
of this system set equal to zero is the equation for natural frequencies. 

In the above procedure, if equation (3.19) Is used Instead of (3.11), wa 


finally obtain the equation for the eigenvalues of the operator L^. 

From equation (3.10) It Is obvious that X^ have the form 


X ■ p (o^ + u ^) 
n n 

and then 

where u are the natural frequencies of thft beam simply supported at 
n 

both ends. The corresponding eigenfunctions are obtained from 
equation (3.16). In summary, a general solution has been presented 
for the free and forced vibrations of a finite Euler- Bernoulli beam. 
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